Abstract. We discuss the Aromatic Infrared Band (AIB) 
profiles observed by ISO-SWS towards a number of bright 
interstellar regions where dense molecular gas is illumi- 
nated by stellar radiation. Our sample spans a broad range 
of excitation conditions (exciting radiation fields with ef- 
fective temperature, T^fj, ranging from 23,000 to 45,000 
K). The SWS spectra are decomposed coherently in our 
sample into Lorentz profiles and a broadband continuum. 
We find that the individual profiles of the main AIBs at 
3.3, 6.2, 8.6 and 11.3 /im are well represented with at 
[ most two lorentzians. The 7.7 ^m-AIB has a more com- 

■ plex shape and requires at least three Lorentz profiles. 
' Furthermore, we show that the positions and widths of 

■ these AIBs are remarkably stable (within a few cm~^) 
O ^' confirming, at higher spectral resolution, the results from 

^ " ISOCAM-CVF and ISOPHOT-S. This spectral decompo- 
\ sition with a small number of Lorentz profiles implicitly 

■ assumes that most of the observed bandwidth arises from 
I I a few, large carriers. Boulanger et al. (1998&) recently pro- 

^ ■ posed that the AIBs are the intrinsic profiles of resonances 
I in small carbon clusters. This interpretation can be tested 
' by comparing the AIB profile parameters (band position 
and width) given in this work to laboratory data on rele- 
vant species when it becomes available. 
Taking advantage of our decomposition, we extract the 
profiles of individual AIBs from the data and compare 
them to a state-of-the-art model of Polycyclic Aromatic 
O ^- Hydrocarbon (PAH) cation emission. In this model, the 
position and width of the AIBs are rather explained by a 
redshift and a broadening of the PAH vibrational bands 
^ ' as the temperature of the molecule increases (Joblin et 
Ci al. 1995). In this context, the present similarity of the 
^ AIB profiles requires that the PAH temperature distribu- 
tion remains roughly the same whatever the radiation field 
hardness. Deriving the temperature distribution of inter- 
stellar PAHs, we show that its hot tail, which controls the 
AIB spectrum, sensitively depends on N^in (the number 
of C-atoms in the smallest PAH) and T^ff. Comparing 
the observed profiles of the individual AIBs to our model 
results, we can match all the AIB profiles (except the 
8.6 /Ltm-AIB profile) if Nmin is increased with T^ff. This 
increase is naturally explained in a picture where small 
PAHs are more efficiently photodissociated in harsher ra- 
diation fields. The observed 8.6 /im-profile, both intensity 
and width, is not explained by our model. 
We then discuss our results in the broader context of ISO 
observations of fainter interstellar regions where PAHs are 
expected to be in neutral form. 
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I. Introduction 

The family of infrared features at 3.3, 6.2, 7.7, 8.6, 

II. 3 and 12.7 /xm has been observed towards a large 
number of sightlines in the Galaxy and in other galaxies 
since the nineteen seventies. Early on, it was recognized 
that these bands correspond to vibrational modes in 
carbonaceous aromatic systems (Duley & Williams 1981, 
Leger & Pugct 1984, Allamandola et al. 1985). These dust 
bands arc therefore dubbed the Aromatic Infrared Bands 
(hereafter AIBs). The bands at 3.3, 8.6, 11.3 and 12.7 
yum stem from vibrational modes of the aromatic C-H 
bond; the remaining bands arc ascribed to vibrations of 
the aromatic C-C bonds. The Infrared Space Observatory 
(ISO, Kessler et al. 1996) mission has provided us with 
an unprecedented wealth of data in this context. 

Comparative studies of the AIBs in a wide variety 
of environments have been made possible by the high 
sensitivity of the ISO camera (ISOCAM, Cesarsky C. 
et al. 1996) with its circular variable filter and of the 
ISOPHOT-S spectrophotometer (Lemke et al. 1996), 
which both have low spectral resolutions (A/AA ^ 40 and 
90 respectively: Boulangcr et al. 1996, 1998a and 19986; 
Cesarsky et al. 1996a, 19966, 2000a and 20005; Crete et 
al. 1999; Klein et al. 1999; Laurcijs et al. 1996; Mattila et 
al. 1996; Pcrsi et al. 1999; Uchida et al. 1998 and 2000). 
The AIB profiles as seen in these data arc very similar 
(in position and width) over a range of objects where 
the stellar radiation field and effective temperature vary 
greatly (1 to 10^ times the standard interstellar radiation 
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field, Te//=11,000 to 50,000 K). The profile invariance 
as well as the large width of the AIBs lead Boulangcr et 
al. (19986) to conclude that the carriers of these bands 
are large aromatic systems containing more than 50 
C-atoms. 

The Short Wavelength Spectrometer (SWS, de 
Graauw et al. 1996) onboard ISO, less sensitive than ISO- 
CAM but with higher spectral resolution (A/AA=1000) 
and broader wavelength coverage (2.4-45 ^m), brings a 
better view of the interstellar AIBs in the brightest re- 
gions. Such detailed data have the ability to constrain the 
nature and physical state of the band carriers (Beintema 
et al. 1996, Molster et al. 1996, Roelfsema et al. 1996, Ver- 
straete et al. 1996, Moutou et al. 2000, van Kerckhoven 
et al. 2000, Hony et al. 2000). In the first part of this pa- 
per, we decompose the AIB spectrum into Lorentz profiles 
and a broadband continuum in order to characterize (band 
position and width) each individual AIB and to compare 
them between objects. In the second part of the paper, we 
confront these new observations with a model considering 
free-flying aromatic molecules (Polycyclic Aromatic Hy- 
drocarbons or PAHs) as the origin of the AIBs. Indeed, the 
presence of AIBs in the low-excitation diffuse interstellar 
medium (Boulangcr et al. 1996, Mattila et al. 1996, Onaka 
et al. 1996) requires the existence of frec-flying PAHs or 
small grains excited by starlight; furthermore, these PAHs 
or grains must be small enough to undergo strong temper- 
ature fluctuations leading to emission of the AIBs in the 
near-infrared (Sellgren 1984). In this emission mechanism, 
the shape of the emergent AIB spectrum only depends on 
the radiation field hardness (or Teff, see Sect. 2) and not 
on the flux of stellar photons (parameterized as x in Sect. 
2). 

In Sect. 2, we present the SWS spectra. The observed 

AIB profllcs are characterized in Sect. 3. These results are 
compared to the predictions of the PAH model in Sect. 4. 
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We summarize and discuss the significance of our results 
in Sect. 5. 



2. Observations and physical conditions 

SWS has observed the AIB spectrum along a number of 
interstellar sightlines covering a wide range of excitation 
conditions. We discuss here three spectra which sample 
well the radiation field sequence covered by the SWS ob- 
servations. The lines of sight selected all correspond to in- 
terface conditions, i. e. , regions where fresh molecular ma- 
terial is directly exposed to the stellar light. At interfaces, 
the AIB emission is usually strong (probably because of 
an enhanced PAH abundance, Bernard et al. 1993): this 
is why we selected such regions to carry out the present 
study. 

At the low excitation end, we have the reflection nebula 

NGC 2023 (TDT=65602309, SWSOl speed 3). This spec- 
trum has already been presented in Moutou et al. (1999). 
SWS looked at a filament, 60" south of the central star 
(central star: RA= 5^ 41"^ 38.3^ DEC= -2°16' 32.6"), 
which is bright in fluorescent H2 emission (Field et al. , 
1998). At the high excitation end, we present here a spec- 
trum of the M17-SW photodissociation interface (RA= 
18*^ 20™ 22.1% DEC= -16°12' 41.3"; TDT=32900866, 
SWSOl speed 4): this is position number 6 of the data 
presented in Verstraete et al. (1996). Finally, wc also have 
the Orion Bar (RA= 5^ 35"^ 20.3% DEC= -5°25' 20"; 
TDT=69501806, SWSOl speed 4) at the position of the 
peak of fluorescent II2 emission (van der Werf et al. 1996). 
All positions given above are in J2000. 

The data reduction was undertaken with the SWS- 
IA3 environment running at the Institut d'Astrophysiquc 
Spatiale, Orsay. The spectrum of NGC 2023 lacks a small 
range around 4 /xm because of bad dark current mea- 
surements. The flux calibration files CAL-G version 030 
have been used. For the beam sizes, we took the values 
recently determined by Salama (2000). This assumes that 
the source completely fills the beam. To check this assump- 
tion, we compared our SWS spectra of M17-SW and of the 
Orion Bar to CAM-CVF data (Cesarsky D. et al. 1996 
and Cesarsky D. et al. 20006): the continuum fiuxes (per 
solid angle) of the two instruments were found to agree 
within 20%. In the case of NGC 2023, the emission seen in 
the ISOCAM-map of Abcrgcl et al. (2000, in preparation) 
looks homogeneous at the position of our SWS spectrum. 
The 6 arcsecond pixels of ISOCAM are much smaller than 
the SWS field of view, and since the ISOCAM image of 
each of our sources is smooth in the region observed by 
SWS, we can safely say that our sources uniformly fill 
the SWS beam. This statement only holds over the 5-16 
fim wavelength range. In fact, to assure continuity in our 
spectra, we had to deviate from the Salama SWS beam 
sizes above 27 /im (the spectral bands 3E and 4 of the 
SWS, see de Graauw et al. 1996). 



Fig. 1. The SWS 2.4-25 /im spectra of the AIBs to- 
wards various interstellar sources (see text). The spectra 
have been normalized to the 6.2 /zm-feature of NGC 2023; 
the scaling factors are given in the figure, e.g., M17-SW/3 
means that the flux of M17-SW has been divided by 3. 
A positive constant has been added to all spectra except 
that of NGC 2023 to ease comparison. The zero flux level 
of each spectrum is indicated by the dashed line. The re- 
solving power (A/AA) is 200 for NGC 2023 and 500 for 
M17-SW and the Orion Bar. The narrow, unresolved lines 
are due to Br/3 2.63 /xm , Bra 4.06 /zm , forbidden ionic 
lines ([Aril] 6.98 /zm, [Arlll] 8.99 /im, [SIV] 10.5 /im, [Nell] 
12.8 /xm, [Nelll] 15.5 /im and [SIII] 18.7/im) and H2 pure 
rotational lines (0-0 S(5) 6.91 /xm, 0-0 S(4) 8.02 /xm, 0-0 
S(3) 9.66 yum, 0-0 S(2) 12.3/im and 0-0 S(l) 17.0 /tm). 

The resulting 2.4-25 /im spectra arc shown in Fig. 1. 
As we will see below, this reduced spectral range amply 
suffices in our discussion of the AIB profiles. The NGC 
2023 spectrum has a resolving power A/AA = 200 and 
those of M17-SW and of the Orion Bar have A/AA = 500. 
The physical conditions prevailing in each interstellar re- 
gion are summarized in Table 1. We assumed that the 
radiation field from the exciting star is a blackbody char- 
acterized by an effective temperature. In this table, we 
give the effective temperature of the exciting star (Teff), 
the dilution factor for the blackbody {Wdu ) , the flux scal- 
ing factor (x) in units of the Habing field (at A — 1000 A) 
and d^, the distance of our line of sight from the central 
star (wc assumed d^, to be the projected distance on the 
sky in all cases). In the case of M17-SW, several stars ex- 
cite the regions we observed: in Table 1 we give effective 
values {T^ff, Ri,) that reproduce well the sum of all radi- 
ation fields. Our excitation sequence thus goes from NGC 
2023 (cool star) to M17-SW (hot star). 

3. The observed AIB profiles 

The continuum below the AIBs shows a clear evolution 
from NGC 2023 to the Orion Bar (see Fig. 1): continuous 
emission of hot, small (radii of a few 10 to 100 A, Desert 
et al. 1990) grains between 10 and 20 /xm is prominent 
in M17-SW and the Orion Bar whereas this is completely 
absent in NGC 2023. This is a consequence of the stronger 
flux and harder radiation field in the Orion Bar which 
heats up these small grains. 

The 6.2 and 7.7 /im bands do not vary much with 
respect to each other whereas they do relative to the 3.3, 
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Table 1. Physical conditions for the AIB spectra - see 
text 



Line of sight NGC 2023 Orion Bar M17-SW 



Teff (10* K) 


2.3" 


3.7= 


4.5" 


R* (R©) 


3 2a,i. 


8.5'' 


12 


d* (pc) 


0.14 


0.24 


1.10 




2.66 


6.38 


0.96 




1.2 


42 


12.5 



The distances to NGC 2023, M17-SW and the Orion Bar were 
taken to be 480 pc", 2.2 kpc*" and 460 pc'' respectively. 
References: (a) Buss et al. 1994, (6) Lang 1991, (c) Rubin et 
al. 1991, (d) van der Werf et al. 1996, (e) FeUi et al. 1984 

8.6 and 11.3 /im AIBs. Such variations cannot be in- 
terpreted as the result of different emission temperature 
distributions of PAHs (obtained with different size distri- 
butions and/or radiation field effective temperature). In 
fact, the 3.3 /um-band is dominated by the emission of 
small (hot) PAHs while the 11.3 /um-band is contributed 
to by larger (cold) PAHs (Schuttc et al. 1993): thus, chang- 
ing the PAH temperature distribution cannot explain why 
the 11.3 /um-band varies along with the 3.3 /zm-band. In- 
stead, modifications of the physical state of PAHs have 
to be invoked. Theoretical (de Frees et al. 1993; Pauzat et 
al. 1995,1997; Langhoff 1996) and laboratory (Szczepanski 
& Vala 1993, Hudgins et al. 1994, Hudgins & AUamandola 
1995, Hudgins & Sandford 1998) studies showed that ion- 
ized species have stronger C-C bands. On the other hand, 
dehydrogcnatcd PAHs have weaker C-H bands (Pauzat et 
al. 1995, 1997). As we show in Sect. 4.4.1, the weaker C-H 
band emission in M17-SW can be explained if the hydro- 
genation fraction is significantly reduced. 

Moreover, the AIB profiles show little spectral sub- 
structure even though most AIBs are fully resolved in 
our SWS data (the resolving power is A/ A A = 200 
to 500). A direct comparison of the spectra displayed 
in Fig. 1 shows that the position and width of the 
major, simple AIBs (at 3.3, 6.2, 8.6 and 11.3 ^m) arc 
roughly the same. The decomposition of the spectra 
we discuss below primarily aims at quantifying this 
comparison and also at disentangling in a systematic way 
the profile of a given AIB from the other bands as well 
as from the underlying continuum due to hot small grains. 

Boulanger et al. (1998&) showed that the AIBs in 
CAM-CVF data can be decomposed into Lorentz pro- 
files and a linear underlying continuum. It must be em- 
phasized, however, that the AIBs are barely resolved in 



the CAM-CVF data (in particular the 6.2, 8.6 and 11.3 
lj,m bands, see Table 2 in Boulanger et al. 1998& and Fig. 
3 of Cesarsky D. et al. 2000). The SWS data at high spec- 
tral resolution presented here do not suffer from this lim- 
itation. The use of a lorentzian band shape implicitly as- 
sumes that the AIB profiles arise from the intrinsic width 
of molecular transitions and/or resonances in small solid 
particles (e. g. , Bohren & Huffman 1983). 

In the current paradigm, the AIBs result of the super- 
position of many vibrational bands produced by a pop- 
ulation of interstellar PAHs with a wide range of sizes 
(molecules containing a few tens to a few hundreds car- 
bon atoms. Desert et al. 1990, Schuttc et al. 1993). In 
this respect, laboratory and theoretical studies on small 
PAH species teach us that the band shapes (position and 
width) of vibrational transitions (i) depend on the tem- 
perature of the molecule (Joblin et al. 1995) and, (ii) vary 
from one species to another (in particular, the position of 
the vibrational bands depends on the size and symme- 
try of the molecule: Szczepanski & Vala 1993, Hudgins & 
AUamandola 1995 and 1998, Joblin et al. 1995, Langhoff 
1996). The observed AIBs may actually result from a com- 
bination of these two effects. Elaborating on these studies, 
variability in the AIBs (changing band ratios, presence of 
substructure in the band profiles) from different interstel- 
lar sightlines is predicted as a consequence of a changing 
PAH population and/or different exciting radiation fields. 
Specifically, the position and width of individual bands are 
expected to vary by a few to several tens of cm~^ from 
one species to another and/or as a result of different emis- 
sion temperatures. Some changes in the AIB profiles have 
been observed towards H ii regions and reflection nebulae 
(Roelfsema et al. 1996, Verstraete et al. 1996, Cesarsky 
et al. 2000a, Uchida et al. 2000, Peelers et al. 2000 in 
preparation). The case of the general (bright) interstellar 
medium is covered below with a quantitative comparison 
of the AIB profiles. 

In this work, rather than establish the "final" AIB pro- 
file parameters, we aim at comparing on the same footing 
the individual AIB profiles under different excitation con- 
ditions. We have therefore decomposed the present spectra 
into Lorentz profiles and a modified blackbody as under- 
lying continuum. We restricted ourselves to the minimum 
number of Lorentz profiles required in order to produce a 
reasonable overall fit and a good representation of every 
individual AIB. 

We used a classical gradient-expansion algorithm with 
analytical partial derivatives and performed the fit in the 
wavenumber space {x = 1/A in cm"'^) over the 2.4-25 
fj,m. wavelength range. In addition to the Lorentz profiles, 
the underlying continuum is fitted simultaneously. For 
the latter, we took a modified blackbody with an emis- 
sivity law proportional to x, the temperature and peak 
brightness of which were the free parameters. The same 
set of fit parameters was used for all objects. Such a set 
was first fixed on the M17-SW spectrum which has a high 
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signal-to-noise ratio and a good feature-to-continuum 
contrast. To fit the 2.4-25 /xm spectrum, twenty Lorentz 
profiles were necessary. Then, the parameter values of the 
Lorentz profiles and of the continuum in M17-SW were 
used as input to fit the other AIB spectra: very good fits 
were obtained by relaxing first the Lorentz amplitudes 
and blackbody continuum parameters, suggesting a 
complete and robust decomposition of the AIB spectrum. 
After adjusting the amplitudes, the profile (position and 
width) and continuum parameters were fine-tuned, simul- 
taneously, over restricted spectral ranges. The centroids, 
widths and amplitudes of the lorentzian fits to the main 
AIBs are given in Table 2. Our fit to the 5-25 /im-AIB 
spectrum is shown in Fig. 2. The blackbody component 
of this decomposition is consistent with the emission 
of warm dust in the mid-infrared and in particular its 
exponential decay (the Wien tail of the blackbody). The 
observed strong variability of the 20 /xm-continuum flux 
(very weak in NGC 2023 whereas strong in M17-SW and 
the Orion Bar) then reflects the varying temperature of 
the warm dust component. As can be seen in Fig. 2, the 
underlying continuum contributes little below the AIBs. 
On the other hand, additional, colder blackbody type 
continua are required to fit the full SWS spectra out to 
45 yum. 

We note that broad bands are required at about 1000 

and 1450 cm"-'^ in order to explain the contimmm be- 
tween the AIBs. These bands may not be associated with 
the AIBs but, for simplicity, we assumed their profiles to 
be lorentzian. Their parameters are not well constrained 
in our decomposition: the sole requirement is that the 
corresponding profiles are broad enough to reproduce the 
smooth continuum observed in these spectral regions. The 
fitted profiles of the neighbouring AIBs are somewhat sen- 
sitive to the widths adopted for the 1000 and 1450 cm~^ 
bands: for instance, if the full width at half maximum 
(FWHM) of the 1450 cm^^ band is increased from 200 to 
300 cm-\ the 6.2 /xm-band has its FWHM reduced by 
2.5 cm""'^ and its position redshifted by 0.8 cm""'^. In or- 
der to coherently compare the AIB profiles, we have fixed 
the width of these broad bands: namely, FWHM = 300 
cm-i for the 1000 cm'^-band and FWHM = 200 cm"! 
for the 1450 cm^^-band. At this stage, we can point out 
that combinations of PAH vibrational modes have been 
predicted to accumulate between 1000 and 2000 cm~^ in 
a broad structure (Bernard et al. 1989). 

Also noteworthy is that two lorentzians are required to 
correctly reproduce the red wing asymmetry of the 1609 
cm~^ (6.2 fim) and 890 cm~^ (11.3 ^m) bands: these 
components are labelled "core" and "red wing" in Table 2. 
The feature centered around 1300 cm^^ (the classical "7.7 
/xm-band") shows 3 sub-peaks at about 1273, 1310 and 
1328 cm~^. In the following, we will call the narrow 1310 
cm"'^ peak the "7.6 /im-band" while the broader 1273 
and 1328 cm~^ components will be dubbed the "7.8" and 



Table 2. Fit parameters of the AIB Lorentz profiles 



AIB NGC 2023 


M17-SW Orion Bar 


11.3 /im 


888.6 1 


889.9 


889.7 




20.8 ^ 


17.8 


14.2 




877 ^ 


2292 


7241 


Core 


889.3 


890.2 


889.9 




13.1 


10.1 


10.9 




564 


1644 


5928 


Red wing 


881.9 


880.3 


880.5 




29.6 


30.4 


30.0 




384 


915 


1820 


8.6 nm 


1164.0 


1163.8 


1161.5 




49.3 


44.6 


48.0 




355 


780 


1957 


7.8 fim 


1274.2 


1273.5 


1275.1 




70.2 


67.6 


54.2 




392 


1494 


2417 


7.6 jim 


1309.9 


1312.6 


1311.3 




22.9 


28.6 


25.6 




130 


912 


1760 


7.5 /xm 


1328.1 


1327.8 


1329.7 




68.2 


74.5 


55.9 




466 


1074 


2022 


6.2 /xm 


1608.1 


1607.7 


1609.7 




48.6 


42.6 


38.9 




550 


1748 


3428 


Core 


1608.5 


1608.6 


1610.5 




17.1 


30.0 


25.5 




224 


1376 


2369 


Red wing 


1600.5 


1586.2 


1594.9 




80.0 


64.4 


64.8 




338 


546 


1291 


3.3 /xm 


3041.8 


3039.1 


3040.0 




43.0 


38.8 


40.4 




98 


171 


834 


center in cm" 


-1 (± 0.8 


cm~^, see text) 




width in cm" 


-1 (± 2 cm-i) 




amplitude in 


MJy/sr 







"7.5 /xm-bands" respectively. The observed 1310 cm~^- 

feature has a narrow core and a broad blue wing which 
demands another component at 1328 cm~^. A lorentzian 
is also required around 1204 cm~^ (FWHM ^ 70 cm~^) 
to fill the gap between the 7.8 and the 8.6 /im-features. In 
the case of Orion this band shifts to 1209 cm^^ in order 
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Fig. 2. The decomposition of the AIB spectra between 
5 and 25 /zm. The resolving power (A/AA) of the data is 
200 for NGC 2023 and 500 for M17-SW and the Orion 
Bar. Each Lorentz component is represented by a dotted 
line, the underlying continuum by a dashed line. The full 
fit is the solid line. 

to reproduce the pronounced and extended blue wing of 
the 8.6 /xm-band. 

In Fig. 3 we show the profile and fit of the 3.3 
/im band. The band shape is also lorcntzian and the 
continuum has the same functional form as in the 5-25 
fim region. In other words, we fitted the 2.4-25 /um con- 
tinuum with a single blackbody component. 

The decomposition presented above is not unique: we 
also tried to use a multicomponent power-law continuum 
{A/x^ + B/x^ + C/x + D with x = 1/A and where 
A, B, C and D are constant parameters) below the AIBs 
from 2.4 to 16 /um. Such a continuum gives a notable con- 
tribution below the AIBs, in particular, it replaces com- 
pletely the 1000 cm~^ broad band. However, a multicom- 
ponent power-law continuum has no straightforward phys- 
ical interpretation and it cannot describe the strong rise 
of the spectrum beyond 16 /um observed in M17-SW and 
the Orion Bar. Nor can it accommodate the weak 20 /xm- 
flux in NGC 2023. In any case, comparing the fits obtained 
with the two types of continua (power-law and blackbody), 
we found the positions and widths of the main AIBs char- 
acterized in Table 2 to vary by less than 1.6 and 4 cm"-'^ 
respectively in a given object. Finally, using gaussian band 
shapes for the AIBs, we got as good a decomposition (see 
also Boulanger et al. 19986): yet the continuum under the 
AIBs was stronger and more structured because of the 
weaker profile wings. We believe that the fit quality does 
not preclude any band profile: in fact, the SWS AIB spec- 
trum is so rich that there are always enough parameters 
in the fit to accomodate any choice of band profile. 

We have thus determined in a coherent way the posi- 
tion and width of the strong, well-delineated AIBs at 3.3, 
6.2, 7.6, 8.6 and 11.3 /xm with an accuracy of 0.8 and 2 
respectively. Inspection of Table 2 shows that there 
are significant variations in the width of the 6.2 and 11.3 
/xm-bands (9.7 and 6.6 cm~^ respectively) as well as in 
the position of the 3.3 and 7.6 /xm-bands (2.7 cm"-'^ in 
both cases). However, all these variations come within the 
accuracy range given above when the values of NGC 2023 
are excluded. In this latter object, the poorer signal-to- 
noise ratio and resolving power (see Fig. 1 and Figs. 8 to 



Fig. 3. Same as Fig. 2 in the region of the 3.3 /xm AIB. 
Note that the continuum has the same functional form 
than that already seen in Fig. 2. In the Orion Bar spec- 
trum, the Pf(5 hydrogen recombination line clearly visible 
on the top of the 3.3 /xm-band was masked to perform the 
fit. The Pf 7 at 3.74 /im and Pf e at 3.04 /im lines can 
also be seen in emission in the Orion Bar spectrum. 

11) has degraded the band profiles: this is probably why 
some AIB parameters are singular in NGC 2023. On the 
other hand, we find that the position of the 8.6 /im-band 
varies by 2.5 cm~^ (Table 2) and that most of this vari- 
ation is due to a redshift in the Orion Bar spectrum. As 
noted above, the 8.6 /im-band in this object has an ex- 
tended blue wing: this spectral change is related to the 
redshift of the band itself. This result may point to the 
more profound modifications seen in this spectral range 
by Roelfsema et al. and Verstraete et al. (1996) towards 
more excited regions. 

From the above discussion, we conclude that most 
AIB profiles (except the 8.6 /xm-AIB in the Orion Bar) do 
not vary to the accuracy of our spectral decomposition. 
At 3.3 /im, a similar result was already obtained by 
Tokunaga et al. (1991) who compared the band profile 
of planetary nebulae and Hll regions (their Type 1 
profile) with excitation conditions (Teff and x) similar 
to that of our sample (see Table 1). Roche et al. (1996) 
confirmed this result on a larger, lower-excitation sample 
of planetary nebulae. Similarly, Witteborn et al. (1989) 
showed that the 11.3 /xm-band profile is rather stable. 
Furthermore, the observed smoothness and invariance of 
the AIBs across a wide range of excitation conditions 
appears difficult to reconcile with the variability expected 
from laboratory results. 

Having performed this mathematical parameteriza- 
tion of the SWS spectra, we are now able to extract 
the individual AIB profiles and to compare them in 
different objects and to model predictions. Yet, what is 
the physical meaning of this spectral decomposition ? We 
have chosen lorentzian band shapes and we represented 
most AIBs (except the 7.7 /xm) with one or two Lorentz 
profiles: this implicitly assumes that the AIBs arise from 
a few vibrational bands common to many PAHs and that 
most of the bandwidth arises from a single carrier. The 
invariance and smoothness of the AIBs is then naturally 
explained (Boulanger et al. 1998&). To verify this inter- 
pretation of the AIBs, the band parameters of Table 2 can 
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be compared directly to laboratory or theoretical studies 
on relevant PAHs (large molecules containing more than 
50 C-atoms as argued by Boulanger et al. 19986). But 
there arc other ways to look at the interstellar AIB 
spectrum which do not result in lorentzian bandshapes. 
For instance, simulations of astronomical spectra based 
on recent laboratory studies of small PAHs (AUamandola 
et al. 1999) have shown that the AIB spectrum may 
be decomposed by assigning a multiplicity of (species 
dependent) vibrational bands to each AIB. Another 
possibility to interpret the AIB spectrum is based on 
the laboratory work of Joblin et al. (1995) which shows 
that the vibrational bands of PAHs are broadened and 
redshifted as the temperature of the molecule is raised. 
In this picture, the profile of a vibrational band from a 
single molecule with a given internal energy has a Lorentz 
shape (the intrinsic profile) which only depends on the 
temperature of the molecule (and not on the species). The 
observed AIBs are then interpreted as the superposition 
of many Lorentz profiles corresponding to all the possible 
temperatures reached by a population of PAHs. 

In the next section, devoted to modelling, we adopt 
and detail this latter view of the AIB spectrum. In this 
context, using the relationship between bandwidth and 
temperature established by Jobhn et al. (1995), we note 
that the width of the observed 3.3 /um-band (40 cm~^) is 
well explained if the emitting PAH has a temperature of 
about 1000 K. Such emission temperatures are in good 
agreement with what is expected for interstellar PAHs 
emitting during temperature fluctuations (sec Sect. 4). Us- 
ing the present spectral decomposition of the SWS data, 
we can now extract the individual AIB profiles and com- 
pare them to the predictions of a PAH emission model 
that uses the best available laboratory data on PAHs. 

4. Probing the PAH model 

To first order, the AIB profiles essentially remain identical 

while the exciting radiation hardens considerably: T^ff 
= 23,000 K in NGC 2023 to 45,000 K in M17-SW corre- 
sponding to the mean energy of the photons absorbed by 
a PAH of 6 and 9 cV respectively. Such contrasted internal 
energies imply differences in the emission temperature of 
a few 100 K for a given molecular size. With this temper- 
ature change, band position shifts and band broadening 
of at least 10 cm~^ are expected (Jobhn et al. 1995). 
How can then the interstellar AIB profiles be so stable ? 
One possibility is that some process keeps the tempera- 
ture distribution of interstellar PAHs essentially the same. 

To investigate this issue, we have computed the tem- 
perature distribution and associated infrared (IR) emis- 
sion of an interstellar PAH population using a modified 
version of the model described in Pcch et al. (2000 here- 
after PJB) where we study in detail the temperature fluc- 



tuations and the temperature distribution of PAHs. Im- 
proving on former work (Leger et al. 19896, Schutte et 
al. 1993, Cook & Saykally 1998), this model takes into 
account the full spectral distribution of the exciting ra- 
diation and implements the recent results of Joblin et 
al. (1995) on the temperature dependence of the band 
proflles of neutral PAHs. In addition, it must be empha- 
sized that the PAH IR emission cross-section of the PJB 
model was derived from studies on PAH cations: indeed, 
we show in Sect. 4.4.1 that singly ionized PAHs reproduce 
better the astronomical AIB spectra confirming previous 
work (Langhofl[ 1996, Cook k Saykally 1998, AUaman- 
dola et al. 1999, Hudgins & AUamandola 1999). This ob- 
servational requirement is supported by theoretical stud- 
ies on the charge state of interstellar PAHs for physical 
conditions comparable to those of our sample of objects 
(Bakes & Tielens 1994, Dartois & d'Hendecourt 1997). We 
first briefly describe the equations giving the PAH cool- 
ing curve and emission, compute the temperature distri- 
bution of interstellar PAHs and then compare our model 
PAH emission spectrum to the data. For the sake of clar- 
ity, we limit the model-to-data comparison to NGC 2023 
and M17-SW which are the extreme cases of our data 
sample and we focus on the well-deflned 3.3, 6.2, 8.6 and 
11.3 /xm-bands. In order to emphasize the temperature 
fluctuations of PAHs and the resulting temperature dis- 
tribution, we have written the model equations in terms 
of cross-sections rather than Einstein A-cocSicicnts: this 
formulation is however strictly equivalent to that of PJB 
(see, e. g. , Schutte et al. 1993, Eq. 15). 

4.1. Emission by PAHs 

In the interstellar medium, a PAH containing Nc carbon 
atoms absorbs radiation from the surrounding stars at a 
rate: 

i?a6s (photons s ) = j j^dv (1) 

with (t" the visible-ultraviolet (UV) absorption cross- 
section of PAHs and F^, the incident stellar flux assumed 
to be a diluted blackbody with the dilution factor and 
eff'ective temperature given in Table 1. The absorption 
cross-section is proportional to Nc', it also has a cut-off 
(due to the electronic transitions) in the visible to near- 
IR range whose wavelength increases with Nc (Verstraete 
& Leger 1992). The mean photon energy absorbed by a 
molecule is given by: 

Eabs = with Pabs = di^- (2) 

J~i-abs JO 

Pabs is the power absorbed per molecule. In the standard 
interstellar radiation field of Mathis et al. (1983), we find 
that a PAH with Nc = 25 or 500 carbon atoms absorbs 
2.4 X 10-27 or 4x10-27 Watt per C-atom with Eabs = 
4.8 or 2.3 eV, respectively. The values for Nc = 25 are 
in good agreement with the laboratory results of Joblin et 
al. (1992) on smah PAHs. 
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As shown by Leger et al. (1989a) and Schutte et 
al. (1993), the IR emission of a PAH is well treated within 
the thermal approximation and using Kirchhoff's law. In 
the following, we will discuss the complete spectrum emit- 
ted by PAHs in terms of the spectral energy distribution 
(SED) noted Si^v) which is equal to vl^. The SED emitted 
by a single molecule containing Nc carbon atoms is thus: 



bs 



where 



v' al, B,, {T,{t,v,Nc)) dtdv (3) 



(4) 



Rabs hv 

is the probability density for absorption of a photon of en- 
ergy hv, all the emission cross-section in the IR, t^. is the 
time at the end of cooling (see Sect. 4.2) and Tc[t, v, Nc) 
is the cooling curve of the molecule. In the following, the 
frequency of photons emitted in the IR will be labelled 
v' and that of stellar photons absorbed in the visible-UV 
range will be written v. 

The PAH IR emission cross-section is defined as in 
PJB. The band profiles have a lorentzian shape and at 3.3, 
6.2, 8.6 and 11.3 ^m, the band position and width follow 
temperature dependent laws. Each Lorentz profile is nor- 
malized to the peak value of the band cross-section, a^ 
and we assume that the integrated cross-section (or Ein- 
stein A-coefficient, see Schutte et al. 1993) of each band 
profile is conserved with temperature. As already noted 
above, the resulting IR cross-section describes the proper- 
ties of singly-ionized PAHs. To match the observed AIB, 
the 8.6 /im-band had to be multiplied by 3. Because the 
available laboratory results do not account for the complex 
shape of the observed 7.7 /um-band (see Sect. 3), we define 
an empirical 7.7 /xm-band with a shape derived from the 
observations (profile centre at 1300 cm"^ and FWHM — 
113 cm""'^ from the present data) and with the (laboratory 
determined) Einstein A-coefficient of PJB. This profile of 
the 7.7 /xm-band is furthermore assumed temperature in- 
dependent. We also added the 16.4 /xm-AIB detected by 
ISO towards many sightlines (Moutou et al. 2000, van Ker- 
ckhoven et al. 2000). This band is often associated with 
the "classical" 3 to 13 ^m-AIBs and its intensity is not 
related to the hot dust continuum appearing in more ex- 
cited objects (see Fig. 1), i. e. , as for the other AIBs. 
We represent the 16.4 /xm-band with a Lorentz profile 
centered at 609 cm^^, with FWHM=6 cm^^ (both be- 
ing observational values, see Moutou et al. 2000) and an 
Einstein ^-coefficient of 3.1 x 10~^ s~^ per carbon atom 
which corresponds to the average value of the laboratory 
measurements given in Moutou et al. (1996, Table 3). As 
we will see below, the intensity of the 16.4 /itm-band can 
constrain the power-law index of the PAH size distribu- 
tion. At longer wavelengths, we adopted the far-IR (A > 
20 /im) cross-section of Schutte et al. (1993). 



The size distribution of PAHs is defined by dN = 
B Nq dNc, i.e., the number of PAHs with a number of 
C- atoms between Nq and Nc + dNc- Writing the size dis- 
tribution with respect to Nc frees us of any assumptions 
on the PAH geometrical shape [|. The number of H-atoms 
per molecule is assumed to be Nh — Jh ^ V^Nc (Omont 
1986, in the case of symmetric PAHs) with /h the hy- 
drogenation fraction of a molecule: if fH=0 the molecule 
is completely dehydrogenated, if fn—^ the molecule is 
fully hydrogenated. Unless otherwise stated, we will as- 
sume fH=^- The mass distribution is then straightfor- 
wardly found as dN x {mc + mn x fHy^/Nc) where mc 
is the carbon atomic mass and mH the atomic mass of hy- 
drogen. The normalization factor, B, is determined from 
the PAH abundance in terms of the fraction of interstel- 
lar carbon ([C/HJ/sm = 2.6 x lO"'', Snow & Witt 1996) 
locked up in PAHs. Finally, the emergent SED emitted by 
PAHs is 

S'{v') = / S{v',Nc)dN{Nc) (5) 

with Njnin and N^^ax the smallest and largest number of 
C-atoms per molecule respectively. We also define the total 
luminosity emitted by PAHs as = /q S^, dlnv' . 



4.2. The infrared cooling of PAHs 

Upon absorption of a UV-photon, the PAH temperature is 
raised to a peak value, Tp. The molecule then cools rapidly 
by the emission of IR vibrational mode photons. The peak 
temperature Tp is found from first principles: 

Uf~U^^ j " C{T) dT (6) 



where C{T) is the heat capacity of PAHs as given by 
Dwek et al. (1998, Eq. A4 of their Appendix), Uij are 
the initial and final molecular internal energies and Ti 
the PAH temperature prior to photon absorption or 
the mean molecular temperature between two photon 
absorptions. We assumed that Uf » Ui with Uf = hv 
the energy of the absorbed photon. Furthermore, we took 
Ti to be the PAH temperature at the end of cooling (see 
the definition below). The heat capacity, as a function 
of temperature, is approximated by a polynomial in T; 
it is also proportional to Nc but independent of Nh, i. 
e. , the contribution of the hydrogen atoms to C(T) is 
assumed to be always the same whatever the size of the 
molecule (Dwek et al. 1998). In this case and for a given 
energy per C-atom in the molecule {hv/Nc), Tp is fixed. 



^ Using the grain radius, a, as a variable the size distribution 
reads dN ^ a" da. We thus have the relationship P = {a — 
d + l)/d, where d is the dimension of the grain: d = 2 for 
planar grains and d = 3 for spherical grains. For instance, 
assuming PAHs are planar (a — 0.9^/ Nc) a value of -3.5 for a 
corresponds to P — —2.25. 
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The coohng curve Tc{t, v, Nq) is obtained by numerical 
integration {vs. temperature) of the energy conservation 
relationship during cooUng, 

Pe (Te) dt = ^C{T,) dT (7) 
where Pe is the total power emitted at all wavelengths by 
the molecule and dT the temperature drop during the time 
dt. This simple differential equation determines a family of 
cooling curves, Tc, as a solution, each cooling curve being 
completely defined by the knowledge of Nc and of the 
initial condition, i. e. , T = Tp at t = 0. 

The power emitted by a PAH at temperature T is: 

/■OO 

Pe{T)= / al, xA-KB,,{T)dv'. (8) 
Jo 

Cooling curves were computed on a grid of 50 points 
with constant temperature intervals. Under these condi- 
tions, we found that the energy conservation requirement: 

/■OO 

/ Pe{Te{t)) dt = hu (9) 

Jo 

is always met within a few percent. 

Each molecule sees its cooling interrupted by the ab- 
sorption of the next photon which is a stochastic pro- 
cess (Desert et al. 1986). Moreover, gas-grain exchanges 
must be taken into account (Rouan et al. 1992, Draine & 
Lazarian 1998) to reliably estimate the temperature of a 
PAH in the low temperature part of the cooling curve. In 
fact, the emission from the cool (but long-lasting) tail of 
the cooling curves will peak in the far-IR/submillimeter, 
far away from the 3 to 20 /im-rangc on which we focus 
here. The detailed treatment of this problem is outside 
the scope of this paper and we adopt the following simple 
criterion to truncate our cooling curves. The mean time 
elapsed between two photon absorptions is 1/ Rats' it is 
for instance 6.9 and 2.6 hours in the radiation field of 
NGC 2023 and M17-SW respectively and for a molecule 
with Nc = 20. Consequently, all cooling curves have been 
truncated at = 1/ Rats to account for the repeated ab- 
sorptions of UV photons by a molecule; te decreases as 
N^^ because Rats ~ Nc- We thus define the lowest tem- 
perature of a PAH heated by a photon of energy hv as 
Ti = Tcite,!^, Nc), i. e. , equal to the temperature at the 
end of the cooling. For instance, in the radiation field of 
NGC 2023 small PAHs [Nc = 20) eventually reach tem- 
perature as low as 15 K at time te whereas bigger PAHs 
{Nc = 500) remain at around 35 K. 

In Fig. 4, we show the cooling curves obtained for 
various energies per C-atom Ec = hv/Nc- As expected, 
the peak temperatures are identical for a given Ec- We 
note that the Tc-curves do not change much for different 
iVc-values while Ec remains fixed: in fact if the emission 
cross-section were only proportional to Nc all the cool- 
ing curves with {hv^Nc) combinations yielding the same 
Ec would merge. Cooling curves with larger Nc have 
higher temperatures because emission in the C-H modes 
is less important (in a fully hydrogenated PAH we have 
Nh/Nc = \/^/Nc). The small variation in the T^-curves 



Fig. 4. Cooling curves for various energy per carbon atom 
(the zero time point has been excluded). The numbers in 
parenthesis give the value of Ec = hv/Nc in meV/C. The 
corresponding {hv, Nc) combinations were as follows: (o) 
£;c=680 meV/C, solid line is (/iz/=13.6 eV, Nc=20), {b) 
453 mcV/C, solid = (13.6, 30) and dot-dash = (9.06, 20), 
(c) 272 meV/C, solid = (13.6, 50) and dot-dash = (5.44, 
20), {d) 136 meV/C, solid = (13.6, 100) and dot-dash = 
(2.72, 20), (e) 27.2 meV/C, solid = (13.6, 500) and dot- 
dash = (0.54, 20) 

at fixed Ec then reflects the decrease of the emission cross- 
section per carbon atom with increasing size. 

4.3. The temperature distribution of interstellar PAHs 

As noted above, the observed similarity of the AIB 
profiles in our sample means that the range of PAH 
emission temperatures does not vary much: to test this 
idea, we determine here the temperature distribution 
of PAHs. Moreover, as we show below, the effect of the 
size distribution and exciting radiation field parameters 
are illustrated in a synthetic fashion in the temperature 
distribution. 

The observed AIB spectrum results principally from 
the superposition of many blackbodies at different tem- 
peratures (Eqs. 3 and 5). Following Eq. 3 it is clear that 
a given blackbody will have a significant contribution 
to the emergent spectrum if its temperature T is high 
and/or if it remains a long time at that temperature: 
a PAH with a temperature between T and T + dT 
(corresponding to the time interval [t,t + dt]) will thus 
contribute to the total emitted luminosity Le (see the end 
of Sect. 4.1) with a weight proportional to Pe{T) x dt or 
the fraction of the total available energy, hv (the energy 
of the absorbed photon), that is dissipated between t 
and t + dt. This definition of the temperature weight 
is different from previous work (Draine & Anderson 
1985, Desert et al. 1986) which only considered the time 
that a given grain spends in the temperature interval 
[T,T + dT], namely dt/dT: the present definition which 
takes into account the emitted power actually reflects the 
contribution of a given temperature to the flnal spectrum. 

Prom Eq. 7, we see that the temperature weight in the 

total luminosity is simply C{T) x dT . Taking into account 
the probability distribution of the exciting photons and 
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the fraction of PAHs containing Nc carbon atoms, we 
define the temperature weight wt as: 

C(Tc)5T ^ , N^+^5lnNc 

hv NpAH 



figS-gif 

Fig. 5. The distribution of PAH temperatures T x p{T) 
normalized to 1 at peak value. In the panels a), b) 
and c), we used the parameters for the radiation field of 
NGC 2023 (see Table 1). The fixed parameters of each 
plot are given in brackets. In a), we vary N^ax and the 
triple-dot dash line shows the effect of assuming that all 
molecules absorb the same photon {Eats = 6.2 eV). Vari- 
ations of Nmin are explored in b). The power-law index, 
(3, is changed in c). In d), we explore the effect of T^ff 
while keeping the dilution factor of NGC 2023 (see text) 
and Nmin = 20. The solid line corresponds to the case 
of the NGC 2023 radiation field. The dot-dash hne in d) 
shows the case of the M17-SW radiation field with Nmin 
= 30, the value required to fit the observed AIB profiles 
(see Sects. 4.4.1 and 4.4.2). 

where Npah the total number of PAHs and 5T is the 
temperature bin in Tc, 5v the frequency bin of the 
exciting radiation field and 5 In Nq the logarithmic bin 
of the size distribution. All these bins have been taken as 
fixed. With this definition we have J2{t,Nc,v) ~ 

To obtain the PAH temperature distribution in a given 
exciting radiation field, wc build the histogram of all cool- 
ing curves for all molecular sizes in temperature bins of 
constant size (AT = 100 K), each temperature T being 
given the weight wt- The density value in the histogram 
is then simply Y^T^fr" wt' ■ The distribution per temper- 
ature interval, p{T) is then found by normalizing the his- 
togram so that p{T)dT = 1. 

From our weight definition, p{T)dT is the fractional 
contribution of all PAHs with temperatures in the range 
[T, T + dT] to the total emitted luminosity L^. 

The PAH temperature distributions, normalized to 1 

at their peak values to facilitate the comparison, arc shown 
in Fig. 5; they all present a maximum between 300 and 
800 K. The shape of these distributions can be understood 
as follows. Assuming that all molecules absorb a photon 
of energy hv, the temperature weights are proportional 
to C(T) X A^^+V At low T, all the C(T)-curves corre- 
sponding to the different TVc's sum up to give wt- As the 
temperature increases, the largest molecules do not con- 
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Fig. 6. The spectral energy densities corresponding to the 
temperature distributions shown in Fig. 5. All the curves 
have been normalized to 1 at the peak value of the 7.7 /xm- 
band (see text). In panels a) to c), the exciting radiation 
field is that of NGC 2023 and we vary the parameters of 
the PAH size distribution with respect to some reference 
values, namely, Nmin =20, Nmax = 200 and /3 = -2.25. 
In d), we look at the effect of changing the effective tem- 
perature of the exciting radiation field while the dilution 
factor is always equal to that of NGC 2023. In each panel, 
the parameters kept constant are given in brackets. 

tribute anymore to wt because they are too cold: this fade 

out more than compensates for the rise of C{T) with the 
temperature and produces a maximum in the temperature 
distribution. 

Fig. 5 illustrates the sensitivity of the PAH temper- 
ature distribution to the size distribution parameters, 
Nmin, Nmax and /3 (Fig. 5a to 5c) and to the hard- 
ness of the exciting radiation field, parameterized as T^ff 
(Fig. 5d). The associated changes in the AIB spectrum 
are shown in Fig. 6. In all these figures (except the M17- 
SW case in Fig. 5d), we used the same dilution factor for 
the radiation field {Wdu = 2.66 x 10"^"^, corresponding 
to NGC 2023): indeed, as long as one is in the regime of 
temperature fluctuations, the spectral distribution of the 
PAH emission spectrum does not depend on the radiation 
field intensity (Sellgren 1984) while its absolute level scales 
with it (see Eq. 3). Moreover, for the purpose of compar- 
ison, all the spectra of Fig. 6 have been normalized to 1 
at the peak value of the 7.7 /zm-band because its profile 
is independent of the PAH temperature (see Sect. 4.1). 

In Fig. 5a, we consider the effect of taking the full 
spectral distribution of exciting photons instead of the 
same mean absorbed photon energy, Eabs (*• e. , Vi, = 
Slhu — Eabs])- In the following, we will refer to the Eabs- 
approximation whenever it is assumed that each molecule 
absorbs the same photon of energy hv = Eabs at the rate 
Rabs in the given radiation field. The i^afes-approximation 
only affects the high temperature tail of the distribution. 
As we will see in Sect. 4.4.2, this tail matters for the 
width of the 3.3 /xm-band. Nevertheless, the rest of the 
AIB spectrum is little affected by the iJafts-approximation. 
Nmax has a weak influence on the temperature distribu- 
tion (Fig. 5b) and its impact, which is most noticeable 
at long wavelengths (A > 10 fim, see Fig. 6a), is easily 
drowned by small changes of /3, the index of the power-law 
size distribution of PAHs (see Figs. 5c and 6c). When /3 
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decreases, the contribution of the smallest PAHs becomes 
more important and the cold tail of the temperature distri- 
bution is somewhat reduced (Fig. 5c): this leads to an en- 
hanced 3.3/11.3 /Ltm-band ratio and to somewhat broader 
bandwidths (Fig. 6c and Figs. 8 to 11). Note that, due 
to this broadening, the band peak value diminishes (see 
the 6.2 /im-band in Fig. 6c) because the band integrated 
cross-section is conserved with temperature. Both Nmin 
and Teff (which determines Eabs, here Tejf = 10,000 to 
45,000 K corresponds to Eats = 4 to 9 eV) affect the hot 
tail of p{T) (Figs. 56 and 5d) which merely reflects the 
hottest part of the cooling curve (Fig. 4). This behaviour 
of p{T) follows from the fact that the PAH peak tempera- 
ture only depends on the energy content per carbon atom 
in a molecule, Eq (see Sect. 4.2) and that a typical value 
for Ec is Eabs /Nmin {Nmin IS close to the mean value of 
Nc because of the steep power- law size distribution). In 
fact, the temperature distribution and the shape of the 
emission spectrum are similar for different {Eabs, Nmin)- 
couples yielding the same energy per carbon atom. As 
expected, an increase of Nmin leads to colder PAHs (Fig. 
56) which contribute more to the long-wavelength bands 
(Fig. 66) while an increase of Teff has the reverse effect 
(see Figs. 5d and 6d). We pointed out earlier that the PAH 
temperature distribution must remain unchanged in order 
to explain the similarity of the observed AIB profiles: we 
show in Fig. 5d that this condition is fulfilled if Nmin is 
raised along with Te//. Specifically, an increase of Te// 
from 23,000 K (the case of NGC 2023) to 45,000 K (the 
case of M17-SW) requires Nmin = 20 to 30, respectively. 
In the next sections, we show that this requirement allows 
to consistently reproduce the overall AIB spectrum as well 
as the individual band profiles. 

4-. 4- Modelling the AIB spectrum 

We now present the PAH IR-emission as computed from 
the formalism described above in the case of NGC 2023 
and M17-SW. We selected these two objects because they 
present the largest contrast in the effective temperature 
of the exciting radiation field: we therefore expect these 
spectra to provide the strongest constraints on our model. 

4.4.1. The complete spectra 

Having determined the temperature distribution of a pop- 
ulation of PAHs, we can now compute the corresponding 
IR emission using eqs. 3 and 5. The PAH size distribu- 
tion parameters were constrained as follows. We note that 
the overall shape of the 3-13 /j,m AIB spectrum which is 
dominated by emission from molecules at T > 400 K is not 
much affected by changes in Nmax or /3 (see Figs. 6a and 
6c). For simplicity, we fix from now on Nmax = 200. First, 
Nmin is fixed so as to reproduce the observed 3.3/11.3 
/Ltm-band ratio in a given exciting radiation field (de- 
scribed by Teff): we require Nmin = 20 in NGC 2023 while 
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Fig. 7. The spectral energy density in Watt per hydro- 
gen for our sample of objects from the SWS data (dotted 
curve) and assuming a column density of 1.8 10^"'^ cm"'^. 
The resolving power (A/AA) is 200 for NGC 2023 and 
500 for M17-SW. The model PAH emission is the dot- 
dash line. The solid line shows the model complete spec- 
trum which includes the broadband continuum described 
in Sect. 3. We took Nmax = 200 and /? = -3.5 for both 
objects. For NGC 2023, we take Nmin ^ 20, fn = 0.8 
and 10% of the interstellar carbon in PAHs while we re- 
quire Nmin — 30, fn = 0.5 and 8% of the interstellar car- 
bon in PAHs for M17-SW (see text). In the upper panel, 
the dashed line shows the case of neutral PAHs from our 

model with Nm^n = 20, Nmax = 200, p = -3.5, fH = 1. 

For clarity, the SED of neutral PAHs has been normalized 
to the peak value of the 11.3 /xm-band in the SED of ion- 
ized PAHs (dot-dash line). This normalization amounts 
to a PAH abundance representing 6% of the interstellar 
carbon. 

Nmin = 30 in M17-SW. The hydrogenation fraction, fn, 
is constrained from the 7.7/11.3 /im-band ratio: a good 
match to the observed band ratios is obtained with fn ~ 
0.8 for NGC 2023 and fn = 0.5 for M17-SW. The index of 
the power law size distribution, /3, can be constrained from 
the requirement that the predicted 16.4 /xm-band (on top 
of the broadband continuum) matches the observed band: 
we find that /3=— 3.5 reproduces the 16.4 /im-band well 
in both objects |^ Unfortunately, individual PAH species 
show a large spread in the integrated cross-section of the 
16.4 /im-band (see Table 3 of Moutou et al. 1996) which 
does not constrain /?: specifically, to match the observed 
16.4 /im-band we require (3 — —2.25 for the lowest cross- 
section value while P — —6 ioi the largest cross-section 
value. More laboratory work on PAHs is warranted in or- 
der to understand the trend of the 16.4 /tm-band profile 
with respect to Nc and to the various isomer states. As 
discussed by PJB, an additional constraint on /3 is pro- 
vided by the intensity and width of the AIB profiles (see 
the next section). The /3- value can then be adjusted in or- 
der to match the peak position and the FWHM of the 3.3 
/im-band (see Sect. 4.4.2): in fact, the value of -3.5 found 
from the 16.4 /xm-band gives a good fit in the case of NGC 
2023 and M17-SW (see Fig. 8). 

^ This result is not sensitive to the amount of cooling oc- 
curring in the far-IR (A < 20 /im): multiplying or dividing our 
far-lR cross-section by 10 did not alter much our match of the 
16.4 /im-band. 
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At this point, it is interesting to note that that the 6- 
to-13 /im part of the spectrum is produced by molecules 
with emission temperatures between 400 and 1000 K. Dur- 
ing temperature fluctuations, most PAHs go through these 
temperatures: this is why the maximum of the tempera- 
ture distribution (Fig. 5) always falls in this range. In fact, 
the emission of any PAH with Nc < 160 shows ratios of 
band peak values which are quite comparable to the ob- 
served 6-13 /um AIB spectrum (see PJB and Draine & 
Li 2001). On the other hand, the match of the observed 
band position and width in the present data sample and in 
CAM-CVF spectra of similar interstellar regions requires 
small molecules (see next section). However, taking into 
account the full size distribution of PAHs will add emission 
mostly at 3.3 /xm (from the smaller species, see Fig. 6b) 
and in the far-infrared (from the larger species, see Fig. 
6a) without changing much the 6-13 /xm region. Thus, 
in the framework of the present model, the observed in- 
variance of the AIB spectrum in CAM-CVF data (6-16 
/xm) implies that the AIB carriers are small enough to 
undergo large temperature fluctuations (at least 400 K) 
during which they emit in this wavelength range. 

Our best fit SED's are displayed on Fig. 7: we 

assumed a total column density of 1.8x10^-'^ cm~^ (i. 
e. , 1 visible magnitude of extinction) for both NGC 
2023 and M17-SW with 10% and 8%, respectively of 
the interstellar carbon in PAH cations (using [C/H]/5m 
= 2.6x10"'', Snow & Witt 1996). The hydrogenation 
fraction of PAHs was taken to be 80% for NGC 2023 and 
50% for M17-SW. To compare in a meaningful fashion 
the predicted PAH emission to the observed spectrum, 
the broadband continuum described in Sect. 3 has been 
added to our model spectrum. This continuum consists 
of the modified blackbody along with the 1000 and 
1450 cm~^ broad profiles (in M17-SW an additional 
broad profile at 600 cm"-'^ is required, see Fig. 2). All 
the main AlBs (at 3.3, 6.2, 7.7, 8.6 and 11.3 /xm) are 
matched within 20%. We note that the model fails to 
reproduce the 12.7 /im-band in M17-SW: this may be 
due to an ill-definition of the underlying broadband 
continuum in this complex spectral region (see Fig. 2). 
The AIB spectrum corresponding to the neutral PAHs 
(we took the IR cross-sections from Leger et al. 19896) is 
shown in the case of NGC 2023 with a PAH abundance 
corresponding to 6% of the interstellar carbon, with 
Nmin = 20, Nmax = 200 and P = -3.5. The iV„i„-value 
was fixed so as to reproduce the position and width of the 
3.3 /im-band (see next section) and it is the same as for 
the cations. In the case of neutral PAHs, the C-C/C-H 
band ratio is reduced by about one order of magnitude 
with respect to the cations. Matching the observed AIB 
with neutral PAHs would require a very low hydrogen 
coverage {Jh < 0.1) which would produce a strongly 
discrepant 8.6 /im-band and also would be at odds with 
theoretical predictions on the hydrogen coverage of PAHs 



(AUain et al. 1996). 

In the following, we compare in detail the profiles of 
the observed 3.3, 6.2, 8.6 and 11.3/im-AIBs to our model 
results: we restrict our comparison to these bands because 
they are well defined in the data and involve a small num- 
ber of Lorentz fit components. 

4.4.2. Individual AIB profiles 

To do a detailed model-to-data comparison, we now 
study the individual AIB profiles. Due to anharmonic 
couplings, the line profiles of PAH vibrational transitions 
are redshifted and broadened as the temperature increases 
and both follow a linear law (Joblin et al. 1995). The 
temperature dependent broadening laws of P.IB used 
in this modelling do not include the rotational width. 
The centroid and width of each band profile are given 
by z/(T) = lyo + XcXT and Aiy{T) = Auo + Xw x T 
respectively. The coefficients of the linear T-laws have 
been measured in the laboratory on a restricted sample 
of small (Nc < 32), neutral PAHs: we assume that the 
temperature dependence of the vibrational band profile 
are the same for all PAH sizes as well as for cations and 
neutrals. To match the position of the observed AIB 
profiles, we had to slightly change the band centers at 
K: namely, for the 3.3, 6.2, 8.6 and 11.3 /^m-bands 
we took vo= 3079 (3076), 1636 (1627), 1169 (1141) 
and 899 (896) cm~^ where the laboratory values are 
in parenthesis. Such shifts are reasonable in view of the 
laboratory results on the spread in band positions for 
different species. 

In order to highlight the match in band width and 
position, our model profiles are compared to the observed 
AIBs in a normalized fashion. For all the results of best 
fit models represented, the peak absolute intensities do 
not deviate by more than 20% from the observations (Fig. 
7). We use the A^mjn -values derived from the complete 
spectra, i. e. , Nmin = 20 and 30 for NGC 2023 and M17- 
SW respectively. Noting that Nmin has a strong effect on 
the PAH emission (Figs. 5c and 6c and Figs. 8c to 11c), 
it is remarkable that the same value allows to match both 
the absolute peak intensity and the profile shape of the 3.3 
/im-band. The value of /3 is then chosen so as to obtain 
a good match of the width of all bands: we find that f3 
= -3.5 is a good compromise. It must be noted, however, 
that the values of N^in and /3 are strongly dependent on 
the adopted profile temperature laws. 

Figs. 8 to 11 show our model representation of the 
AIB profiles. The contribution of vibrational hot bands 
are not included in our model profiles: in the case of the 
6.2 and 11.3 /im- bands, hot bands provide an additional 
broadening of less than 5 cm~^ (PJB). On Fig. 8c, we 
note the effect of making the i<^af,s-approximation [i. e. , 
all molecules are heated by the same energy photons, see 
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Fig. 8. The normalized profile of the 3.3 /um-AIB as seen 
by ISO-SWS compared to our model emission profile also 
normalized to one (see text). Note that the 3.4 /xm-band 
as well as the underlying plateau (as represented in Fig. 3) 
have been subtracted from the 3.3 /xm-AIB. The resolving 
power of the data (A/AA) is 200 for NGC 2023 and 500 for 
M17-SW. For all the model profiles, we took N^^x = 200. 
In a), the SWS spectrum (from which the continuum and 
the; contribution of all bands other than the 3.3 /im have 
been subtracted) is the dotted curve, the result of our 
model best fit with Nmin = 20 and /3 = —3.5 is the solid 
line and the dashed line represents the ease where f3 = 
—2.25 (while Nmin = 20). In b), we show the contribution 
of the small (20 < Nc < 35, dashed line) and large (35 < 
Nc < 200, dot-dashed Hne) PAHs to the model best fit 
of NGC 2023. In c), the observed 3.3 /Ltm-band profile of 
M17-SW (dotted curve) is shown as well as our model best 
fit (solid line) with Nmin = 30 and /3 = —3.5. The model 
profile for Nmin = 20 (while /3 = —3.5) is shown as the 
dashed line. The dot-dashed line shows a model using the 
-E'obs-approximation {Eats = 9.1 eV, see Sect. 4.3) with 
Nmir-, = 30 and 3 = -3,5. 



Fig. 9. Same as Fig. 8 for the 6.2 /im-band. In the case of 
M17-SW, the profile for the E^abs-approximation has not 
been represented. 

Sect. 4.3): the lack of the hot tail in the PAH temper- 
ature distribution yields a narrower band profile. Hence, 
using this approximation would lead to an underestimate 
of Nmin- e.- .9- . in the case of NGC 2023, Nmin = 15 
is required to match the observed 3.3 /um-band if one 
makes the Sobg-approximation. This approximation has, 
however, negligible consequences for the other bands. The 
effect of a larger /?, -2.25 (corresponding to the classical 
power-law exponent of -3.5 when the size distribution is 
expressed in terms of the grain radius, Mathis et al. 1977) 
is shown as the dashed line on Figs. 8a to 11a: a low /?- 
value is clearly required to match the observed band pro- 
files. Decreasing /3, however, has a moderate impact on 
the band profiles because the then enhanced contribution 



I figlO-Kifl 

Fig. 10. Model results for the 8.6 /xm-band. The data is 
displayed as in Fig. 8 with the following changes. In a), 
the solid line shows the model profile if the bandwidth 
temperature dependence derived from laboratory results 
is used. The dot-dash line shows the model profile for a 
hypothetical broadening 5.5 times faster (see text) and 
the dashed line corresponds to (3 = —2.25 case for the fast 
broadening law. In b), we show the contribution of small 
and large PAHs in the case of the fast broadening law. 
Finally, the display in c) is the same as in a for the M17- 
SW case except for the dash line where Nmin = 20 and 
(3 = —3.5 with the fast broadening law. The unresolved 
lines at 1035, 1110 and 1245 cm"! correspond to 0-0 S(3) 
of H2, a forbidden hne of [Arlll] and 0-0 S(2) of Hg. 



I figll-Kifl 

Fig. 11. Same as Fig. 9 for the 11.3 /xm-band. Note that 
the 11 /um-AIB has not been withdrawn from the data. 

of small molecules is compensated by a strong increase of 
the band width leading to a decrease in intensity (the inte- 
grated cross-sections of the IR bands, a x Au, are assumed 
to be conserved with temperature). 

The harder radiation field of M17-SW results in PAHs 
at higher temperatures (see Fig. 5d), consequently the 
model AIB-profiles are clearly too broad and redshifted 
(see the dashed line in the bottom panel of Figs. 8 to 11) 
if we use the same size distribution as for NGC 2023. For 
the profile width and position, Nmin = 30 is also required 
in order to match the observed AIBs of M17-SW. This 
change of Nmin has a physical interpretation: the photo- 
dissociation and fragmentation rates rise steeply with the 
molecular temperature (Leger et al. 1989a). The Nmin val- 
ues we find are smaller than the theoretical predictions of 
AUain et al. (1996). 

Except for the 8.6 /xm-band, all AIB profiles are well 
explained by the PAH emission model. The hot (small) 
molecules dominate the emission profiles, in particular in 
the AIB wings; they also produce the observed red asym- 
metry (6.2 and 11.3 yum-bands). The cool (large) molecules 
contribute the blue wing of the AIB. We want to stress 
that this match of the AIBs is obtained by using the cross- 
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sections and profile temperature dependence measured on 
small molecules: no additional broadening, reflecting the 
decrease of the level lifetime with increasing molecular size 
(see Sect. 5), has been taken into account here. 

The observed width of the 8.6 /xm-band is not at all re- 
produced by the temperature broadening law derived from 
laboratory work. In fact, a broadening 5.5 times faster 
{Xw = 6.0 10~^ cm~^/K) with temperature would match 
the observations (see Fig. 10). In this latter case, we note 
that (i) the symmetrical shape of the 8.6 /im-band is well 
explained by the slow rcdshift of the band centroid with 
temperature {xc = 0.84x10"^ cm~^/K whereas the other 
AIBs have 2 to 3x10^^ cm^^/K) and (ii) the emission of 
large molecules dominates the core of the 8.6 /xm-band 
(conversely to the other AIBs). 

We have not modelled the 7.7 /im-band because it in- 
cludes several components whose assignment to features in 
laboratory spectra of PAHs is not straightforward. How- 
ever, it is interesting to note that the observed width (23 
to 29 cm~^) of the 7.6 /xm-band (see Sect. 3) is compa- 
rable to the width predicted by PJB. 

We emphasize here that the model profile of a PAH 
emission band is the result of the superposition of many 
Lorentz profiles conversely to the implicit assumption 
behind the decomposition performed in Sect. 3. Consider 
for instance a PAH containing Nc carbon atoms which 
has been heated by a photon of energy hv: the emerging 
band profile of this molecule will be the sum of all the 
Lorentz profiles corresponding to the temperatures of 
the cooling curve, Tc{t, v, Nc) {i- e. , the time integral 
in Eq. 3). The total band profile of a PAH interstellar 
population is eventually obtained from the weighted sum 
of the emerging band profiles from all molecules, for all 
energies of the exciting photon (see Eqs. 3 and 5). 



5. Discussion and Summary 

We have decomposed the ISO-SWS spectra of a number 

of objects covering a wide range of excitation conditions 
in order to study the individual profiles of the Aromatic 
Infrared Bands (AIBs) between 3 and 13 /um. All spectra 
have been decomposed coherently into Lorentz profiles 
and a broadband continuum. We find that the individual 
profiles of the main AIBs at 3.3, 6.2, 8.6 and 11.3 /xm are 
well represented with at most two lorentzians. The 7.7 
/xm-AIB has a more complex shape and requires at least 
three Lorentz profiles. We find that the AIB positions 
and widths are stable to within a few cm"-'^ (see Table 
2) over a range of radiation field hardness (Te//=23,000 
to 45,000 K) thus confirming results gathered with data 
at lower spectral resolution (Boulangcr et al. 1998a, 
Uchida et al. 2000, Chan et al. 2000). This spectral 
decomposition with a small number of Lorentz profiles 
implicitly assumes that, (i) the AIBs arise from a few 
vibrational bands common to many carriers and, (ii) 



that most of the bandwidth arises from a single carrier. 
Boulanger et al. (19986) recently proposed that the 
AIBs are the intrinsic profiles of resonances in small 
carbon clusters containing more than 50 C-atoms. This 
interpretation can be tested by comparing the AIB profile 
parameters (band position and width) given in this work 
to laboratory data on relevant species when it becomes 
available. 

This spectral decomposition, which separates the AIB 
emission from the underlying continuum, allowed us to 
do a detailed comparison of the observed AIB spectrum 
with the predictions of the PAH model where the AIB 
carriers are free-flying aromatic molecules emitting during 
temperature fluctuations. This model uses recent labora- 
tory data and assumes that PAHs are predominantly in 
cationic form as is expected from theoretical work (Bakes 
& Tielens 1994, Dartois & d'Hendecourt 1997). Within 
this framework, the position and width of the AIBs are 
rather explained by a redshift and a broadening of the 
PAH vibrational bands as the temperature of the molecule 
increases (Joblin et al. 1995). The observed similarity in 
the AIB profiles thus requires that some process renders 
the temperature distribution of PAHs rather constant in 
the interstellar regions considered here. Wc first derived 
the temperature distribution for a population of inter- 
stellar PAHs. In particular, we show that the hot tail of 
the tcmperatiire distribution of PAHs (which determines 
the AIB spectrum) depends sensitively on Nmin and T^f / 
which are respectively the size of the smallest PAH (in 
terms of the number of C-atoms in the molecule) and the 
effective temperature of the exciting radiation field. The 
size of the largest PAH, N^ax (by number of C-atoms), 
and the index of the power law size distribution (expressed 
in terms of the number of C-atoms per molecule), /3, were 
found to have little impact on the overall AIB emission. 

We compared our model results to the data in the two 
extreme cases of our sample: NGC 2023 (Te// =23,000 
K) and M17-SW (Te//=45,000 K). We are able to repro- 
duce the spectral distribution of the AIB emission with 
N^,n = 20 in NGC 2023 and iV„i„ = 30 in M17-SW and 
a PAH abundance amounting to 10 and 8% of the inter- 
stellar carbon (these latter values arc in good agreement 
with the constraints set by the extinction curve, Joblin et 
al. 1992, Verstraete et al. 1992, and the infrared emission. 
Desert et al. 1990 and Dwek et al. 1998). The minimum 
PAH size Nmin was found from the requirement that the 
observed 3.3/11.3 /xm-band be well matched. This change 
in Nmin may reflect the enhanced photodcstruction rate 
of small PAHs in regions with harder radiation fields 
(small PAHs reach higher temperatures and evaporate 
efficiently). Using the same A^„ij„-valucs, we find that 
all AIB profile shapes, except for the 8.6 /um-band, can 
be explained with the temperature dependence of the 
band position and width measured in the laboratory. We 
also show that the PAH size distribution must be very 
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steep (/? = -3.5) in order to account for the observed 
bandwidths: in our model representation, the AIB profiles 
are thus mostly contributed to by small PAH species. 
The case of the 8.6 /im-band is anomalous: first, the 
laboratory cross-section is 3 times too weak to explain 
the observed AIB ratios; in addition, we find that the 
broadening of this band with the temperature should be 
5.5 times faster in order to match the profile width seen 
in the SWS data. More laboratory work is required to 
explain this issue. 

In summary, using the present best knowledge of 

PAH spectroscopy, it is possible to account for the AIB 
spectrum towards bright interstellar objects where the 
fraction of singly-ionized PAHs is high. We emphasize 
that both the spectral shape of the AIB spectrum as well 
as the individual AIB shapes are explained consistently 
with a single set of parameter {Nmin and f3) values. In 
this context, the remarkable stability of the AIB profiles 
arises if the hot tail of the PAH temperature distribution 
remains essentially the same, whatever the exciting ra- 
diation field (the low-temperature side corresponding to 
the larger molecules is practically unaffected by changes 
in the radiation field). This requirement is naturally met 
while considering the photodestruction of PAHs: when 
the energy per bond (which is directly proportional to 
the molecular temperature) is sufficient the molecule 
efficiently loses its atoms. A consequence of this is that 
interstellar PAHs would predominantly be destroyed 
by thermal evaporation (photo-thermo-dissociation, see 
Leger et al. 1989a) rather than non-equilibrium processes 
like direct photodissociation (Buch 1989) or Coulomb 
explosion (Leach 1989). 

These results must however be placed in the broader 
context of ISO observations which include less irradi- 
ated regions (x = 1 to 1000). Indeed, ISOCAM-CVF and 
ISOPHOT-S spectra of faint AIB emission highlight the 
pre-requisites of this work. Namely, two strong assump- 
tions have been made in the present model: 

1. the behaviour of the band profiles with the tempera- 
ture (measured on small molecules only, Nc < 50) is 
the same for all PAH sizes. This assumption applies 
for all the AIBs except the 3.3 : in that case the 
band profile is actually dominated by the contribution 
of small species (see Fig. 86). 

2. All PAHs have the infrared emissivity of singly ionized 
cations. 

Assumption 1. is questionable because the intrinsic 
bandwidth (Avq, see Sect. 4.4.2) depends on the density 
of states and on the rate of internal conversion processes 
(e. g. , Leger et al. 1989a). In particular, at a given 
internal energy of the molecule, the density of vibra- 
tional states is known to increase steeply for larger and 
larger molecules (e. g. , Cook & Saykally 1998). Many 



isoenergetic molecular levels may then couple together 
and eventually lead to a decrease of the level lifetime 
(Smalley 1983). Hence, relaxing the first assumption will 
probably result in a significant broadening of the model 
band profiles; this will require a larger minimum PAH 
size in order to reproduce the observed width of AIBs. 
Increasing Nmin, the 3.3 /im-band emission will decrease 
rapidly because larger (colder) molecules emit more at 
longer wavelengths and we will then fail to match the 
intensity of the observed 3.3 /im-band. This problem will 
not be alleviated by assuming a PAH ionized fraction of 
less than 1, [i. e. , not all PAHs are cations): neutral 
PAHs, which have strong C-H/C-C band ratios, (typically 
10 times larger than in cations, Langhoff 1996, Pauzat et 
al. 1997, Hudgins & Sandford 1998), will not reproduce 
the observed 11.3/7.7 /im-band ratio. This relates to the 
second assumption of the present model. 

We have assumed that PAHs are all singly ionized 
(cations) throughout our data sample. Such a choice is 
required to reproduce the observed band ratios, in par- 
ticular the C-H/C-C ratio (Langhoff 1996, AUamandola 
et al. 1999). Yet, this requirement is not fully consistent 
with theoretical predictions (Dartois et al. 1997) of the 
PAH ionized fraction along our lines of sight, namely, 0.5 
(50% of PAHs are cations). Furthermore, this problem 
becomes even more acute in the context of other obser- 
vations spanning a broader range of physical conditions: 
ISOCAM-CVF and ISOPHOT-S data (Boulanger et 
al. 1998a and 1999; Chan et al. 2000; Miville-Deschenes 
et al. 1999; Onaka et al. 1999 and 2000, Onaka 2000; 
Uchida et al. 1998 and 2000) obtained towards regions 
with X = 1 to 10^ (from the diffuse interstellar medium to 
Hll region interfaces and reflection nebulae) show that 
the C-H/C-C ratio of the AIBs is roughly constant and 
always corresponds to that of PAH cations. Over such 
a large range of UV radiation flux, the state of PAHs is 
actually expected to change from fully neutral to fully 
ionized. The ionized fraction of PAHs is determined by 
the value of 7 = x \/T /ue with T the gas temperature and 
Ue the electron density (this parameter is proportional 
to the ratio of the ionization rate to the recombination 
rate of PAHs, Bakes & Tielens 1994). To keep the AIB 
band ratios constant, and hence the ionized fraction of 
PAHs, requires that the ratio ^/T/rie can vary over four 
orders of magnitude in order to compensate the variation 
of X- In reality, T probably does not change by more 
than 2 orders of magnitudes: then rig, which reflects the 
density variations mostly, would have to vary by three 
orders of magnitude across the regions observed to keep 
7 constant. Such large density contrasts at large scale are 
at odds with what is currently known of the structure of 
these interstellar regions. 

More generally, the predicted strong variability of the 
physico-chemical state (ionization, dehydrogenation) of 
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PAHs in space is not reflected in the recent ISO spec- 
troscopy database which covers a variety of astrophysical 
conditions. This prediction was based on studies of small 
PAHs {Nc < 50), the only species currently accessible 
for investigation either in the laboratory or with quan- 
tum chemistry. As shown by Schutte et al. (1993), the 6 
to 16 /.tm-AIB spectrum is also contributed to by large 
PAHs {Nc > 100) while the 3.3 /zm-AIB is dominated 
by the smallest species {Nc < 35, see Fig. 8): this other 
solution (large "PAHs" or carbon clusters) to the 6 to 16 
/xm-AIB spectrum may alleviate the present difficulties. 
However, the spectroscopical and structural properties of 
carbon clusters are poorly known at present: the ISO spec- 
troscopic database is ideally suited to identify plausible 
candidates and stimulate future work on the physics and 
chemistry of carbon clusters. 
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